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Abstract 

Core-valence basis sets for the alkali and alkaline earth metals Li, Be, Na, Mg, K, and Ca are 
proposed. The basis sets are validated by calculating spectroscopic constants of a variety of di- 
atomic molecules involving these elements. Neglect of (3s, 3p) correlation in K and Ca compounds 
will lead to erratic results at best, and chemically nonsensical ones if chalcogens or halogens are 
present. The addition of low-exponent p functions to the K and Ca basis sets is essential for 
smooth convergence of molecular properties. Inclusion of inner-shell correlation is important for 
accurate spectroscopic constants and binding energies of all the compounds. In basis set extrapo- 
lation/convergence calculations, the explicit inclusion of alkali and alkaline earth metal subvalence 
correlation at all steps is essential for K and Ca, strongly recommended for Na, and optional for 
Li and Mg, while in Be compounds, an additive treatment in a separate 'core correlation' step 
is probably sufficient. Consideration of (Is) inner-shell correlation energy in first-row elements 
requires inclusion of (2s, 2p) 'deep core' correlation energy in K and Ca for consistency. The latter 
requires special CCVnZ 'deep core correlation' basis sets. For compounds involving Ca bound to 
electronegative elements, additional d functions in the basis set are strongly recommended. For 
optimal basis set convergence in such cases, we suggest the sequence CV(D+3d)Z, CV(T+2d)Z, 
CV(Q+g?)Z, and CV5Z on calcium. 
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I. INTRODUCTION 



The rate-determining step in accurate wavefunction ab initio calculations is the determi- 
nation of the correlation energy. Moderately reliable methods like CCSD (coupled cluster 
with all single and double substitutions Q]) have theoretical CPU time scalings oc n 2 N 4 , 
where n and N represent the numbers of electrons and basis functions, respectively. For the 
CCSD(T) (i.e^ CCSD with a quasiperturbative estimate of the effect of connected triple 
substitutions |2j) method — which has been termed the "gold standard of quantum chem- 
istry" (T. H. Dunning) — the corresponding scaling is oc n 3 iV 4 . 

For this reason, the idea of only including the 'chemically relevant' valence electrons in the 
correlation problem and constraining the inner-shell or 'core' orbitals to be doubly occupied 
gained currency early on. Bauschlicher et al.|{| were the first to propose a partitioning of 
the inner-shell correlation energy into a core-core (CC) component — involving excitations 
exclusively out of the inner-shell orbitals — and a core- valence (CV) component — involving 
simultaneous excitations out of the valence and inner-shell orbitals. While, at least in 
elements with a large core- valence gap (see below), the CC component can be expected to 
roughly cancel between a molecule and its separated atoms, the CV component may affect 
chemically relevant molecular properties to some degree. 

Over time, it has become recognized that the inclusion of inner-shell correlation — even 
in first-row system - is im po rt an t for .curate biudiug energies (e.g. Q Q), molecular 
geometries (e.g.jj]), and harmonic frequencies (e.g. jjj). The treatment of inner-shell cor- 
relation requires special basis sets which not only have additional radial flexibility, but also 
include 'hard' or 'tight' (i.e., high-exponent) d,f,... functions in order to cover angular 
correlation from the inner-shell orbitals. Special basis sets of this type have been developed, 
such as the cc-pCVnZ (correlation consistent polarized core-valence n-tuple zeta) basis sets 

n fi n 

of Woon and Dunningp for B-Ne, the MT (Martin- Taylor |7() and MTsmall 8] basis sets 
for Li-Ar, the rather small G31arge basis set|2( of Pople and coworkers for the main group 
elements of the first three rows, and very recently the cc-pwCVnZ (i.e. core- valence weighted 
cc-pCVnZ) basis sets of Peterson and Dunning |lfj| for B-Ne and Al-Ar. 



While the CC+CV contribution to absolute correlation energies may rival or exceed the 
valence contribution, its differential contribution to the molecular binding energy is gener- 
ally small compared to the SCF and valence correlation contributions (typically less than 
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1% for lst-row systems). For that reason — as well as the formidable cost of core corre- 
lation calculations — computational thermochemistry schemes (e.g. G3 theory [9|, W1/W2 
theory j^, ^|) that account for inner-shell correlation generally treat the latter as a small 
additive contribution obtained with a relatively compact core correlation basis set. (The 
finding [^] that connected triple substitutions are surprisingly important for CV contribu- 
tions ensures that the core correlation step often is the rate-determining one for benchmark 
thermochemistry calculations, particularly beyond the first row of the Periodic Table.) Re- 
cently, some attempts have been made to replace the core correlation step by bond additivity 
approximations [12, and references therein], and promising results have been obtained 3| by 
means of core polarization potentials [l^j]. 

Despite repeated warnings in the literature against the practice (e.g., by Taylor by one 
of the present authors ^J, and by Woon and Dunning [6]), core correlation is often included 
- for technical reasons or 'because it cannot hurt' — in calculations using basis sets (e.g. 
the standard D UIming ee rrel at 10 n consist ba S1S MU, which are o f mUUm a, b a S1S set 
quality in the inner-shell orbitals) which are not adapted for nonvalence correlation. This 
can often cause errors on computed properties well in excess of the basis set incompleteness 
error: For instance, a comparison between valence-only and all-electron CCSD(T)/cc- 
pVTZ harmonic frequency calculations on the cyclic isomer of C4 suggests core correlation 
contributions to the harmonic frequencies of up to 50 cm -1 , a factor of five more than the 
true correction obtained jlq] with a core correlation basis set. 

For the elements B-Ne and Al-Ar, the gap between valence and inner-shell orbital energies 
is large enough that a conventional core-valence separation of the orbital energies is easily 



made. Things are less simple further down the periodic table. Bauschlicher|20j, |2l| noted 
that for gallium and indium halides, the valence (2s) orbitals of the halogen are below the 
Ga(3d) or In(4d) 'core' orbitals in energy. At the very least, these d orbitals should be added 
to the correlation space; suitable basis sets for this purpose have been developed by Martin 



and Sundermann 



221 ] for Ga and Ge, and by Bauschlicher 



2l| for In. 

But similar issues arise with Group 1 and 2 elements. Radom and coworkers^. Q| 
recently noted catastrophic failures of G2 theory and related methods in predicting the 
heats of formation of various alkali and alkaline earth metal oxides and hydroxides: for 
example, the binding energy of K 2 is underestimated by no less than 256 kJ/mol! Upon 
correlating all electrons, this error decreases to 60 kJ/mol, which is further decreased to 17 
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kJ/mol when more sophisticated electron correlation methods are used. Inspection of an 
atomic orbital energy table quickly reveals why: the oxygen 2s and 2p valence orbitals 
lie below the potassium 3s and 3p inner-shell orbitals, rendering the conventional separation 
between 'valence' and 'core' orbitals essentially meaningless. Inclusion of Group 1 and 2 
subvalence correlation — termed 'relaxed inner valence' (RIV) by Radom and coworkers j3] 
- may in fact be the appropriate treatment, not just for K (and Ca) compounds, but for 
alkali and alkaline earth metal compounds in general. (One of the first papers to recognize 
the importance of inner-shell correlation for alkali metals may have been that by Liu and 
coworkers. 

This problem is far from academic, considering the great importance of K + , Ca 2+ , Na + , 
and Mg 2+ complexes in molecular biology. In addition, a number of these interactions (e.g., 
of the cation- 7r type) are not necessarily amenable to density functional treatments without 
validation of some kind by means of ab initio methods. Clearly, the availability of high- 
quality core-valence basis sets for the elements Li, Be, Na, Mg, K, and Ca would benefit 
the high-accuracy computational thermochemistry and spectroscopy communities as well as 
biomolecular modelers. The purpose of the present work is the development and validation 
of CVnZ (core-valence n-tuple zeta, n=D, T, Q, 5) basis sets for these elements. 

II. COMPUTATIONAL DETAILS 

All electronic structure calculations carried out in the present work were carried out using 



the MOLPRO 2000.1 and 2002.3 program systems [23, |2jj running on SGI Origin 2000 (12 
xMIPS R10000, 300 MHz and 4 xMIPS R10000, 195 MHz, IRIX 6.5) and Compaq ES40 
(4 x EV67, 667 MHz, Tru64 Unix 4. Of) minisupercomputers, as well as Compaq XP1000 
(EV6, 500 MHz, Tru64 Unix 4.0f) and dual Intel Xeon (1.7, 2.0 and 2.4 GHz, RedHat Linux 
7.2 and 7.3) workstations. 

Unless explicitly noted otherwise, energy calculations were carried out at the CCSD(T) 
level; for open-shell systems, single-determinant ROHF reference functions were used, but 
the definition of the CCSD(T) energy according to Watts et al. [29( was employed throughout. 

Validation calculations on diatomic molecules were carried out as follows. Energies were 
computed at 21 equidistant points around the putative bond distance with interpoint spac- 
ings of 0.01-0.03 A, depending on the curvature of the surface. Energies at these points 
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were converged as precisely as feasible, with integral evaluation thresholds being tightened 
as necessary. Polynomials of increasing order were then fit through the points, and the 
significance of the additional order terms subjected to a Fisher-Snedecor test. The expan- 
sion was truncated at the point where the two-way significance of the next higher order 
dropped below 99%. The retained expansion was generally of order between six and eight. 
The minimum of this curve was then sought by means of a Newton-Raphson method, and 
the curve re-expanded around the minimum. A Dunham analysis [3] was then carried out 
on the final curve. The computed dissociation energies at K, D , reported in the tables 
include anharmonic zero-point corrections from the computed uo e and uj e x e at the same level 
of theory. 

As for the accompanying basis sets, standard cc-pVnZ basis sets were used throughout 
on H. For O and F, aug-cc-pVnZ basis sets[31] were used in valence calculations and aug-cc- 
pCVnZ basis setsjy] in calculation where the O and F (Is) cores were correlated. For S and 
CI, the aug-cc-pV(n + d)Z basis sets of Wilson et al.[32j were used in valence calculations, 
and the cc-pCVnZ basis sets of Peterson in calculations where the (2s, 2p) cores of these 
elements were correlated. 

Where deemed necessary, scalar relativistic effects were assessed by means of CCSD(T) 
calculations within the Douglas-Kroll-Hess approximation (33), [mJ as implemented in MOL- 



PRO 2002.3 



281. 



III. OPTIMIZATION OF CVnZ BASIS SETS 



A. General procedure 

The basis set optimizations were carried out using a Fortran program developed in-house. 
Derivatives were obtained numerically using central differences, yielding both the gradient 
and the diagonal of the Hessian. The step thus obtained — essentially Newton-Raphson 
with neglect of off-diagonal Hessian elements — was combined with a Brent line search. Our 
experience with this simple but fairly robust optimization algorithm is quite good as long as 
optimization parameters are not too strongly coupled. In the presence of significant coupling, 
DIIS (Direct Inversion of the Iterative Subspace) [3] dramatically sped up convergence. 

Finite difference step sizes were made roughly proportional to the parameters themselves, 



5 



and increased as necessary if the parameter surface in the affected direction was found to 
be particularly 'flat'. 

Some additional optimizations were carried out using an adaptation of the DOMIN 
program by P. SpellucciQ], which is an implementation of the BFGS (Broyden-Fletcher- 
Goldfarb-Shanno) variable-metric method. Numerical derivatives of order two, four, and 
six were used: the lower orders until an approximate minimum was reached, after which 
the optimization was refined using the higher orders. In most cases, however, the simple 
procedure outlined above appeared to be more robust. 

For the purposes of the optimization (and since only very small systems are involved), 
integral evaluation, SCF and CI convergence thresholds in MOLPRO were tightened to 
essentially machine precision. Basis set parameters were converged to at least four significant 
digits, and five where at all possible. 

In cases where many primitives of a particular angular symmetry are required, we initially 
constrained these functions to follow an even-tempered sequence Q = a^ % ~ n ~ 1 ^' 2 , where n 
is the total number of primitives. Only when optimum a and j3 had been reached were the 
individual exponents optimized further without any constraints. If a set of (n — 1) primitives 
of that symmetry was already available (in casu, from a previously optimized smaller basis 
set), the geometric mean of their exponents was taken as the starting value for a, and the 
arithmetic mean of the ratios between successive primitives as the initial value for (3. 

A few 'checks and balances' can generally be applied to the final exponents: 

1. within a row, the geometric mean of the exponents should increase roughly as (Z*) 2 , 
where Z* is the 'shielded' nuclear charge (first Ahlrichs- Taylor rule f 3 71]): 

2. exponents in a larger set (e.g. 3d) should mesh with the next smaller set (e.g. 2d); 

3. the ratio of successive exponents in a set should be about two or higher, and all gaps 
within the core or valence parts should roughly be of the same order; 

4. the geometric mean of the exponents for each angular momentum should roughly 
increase by a factor of 1.2 for each step up in L (second Ahlrichs- Taylor rulejszl). 
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B. Li and Be 



For Li and Be, the 'unofficial' cc-pVnZ basis sets of Feller [38] supplied in the MOLPRO 
2000.1 basis set library were taken as starting points. The quantity optimized for is the (Is) 
correlation energy at the CISD level for the Li( 2 S') and Be( 1 S') ground states, respectively. 
This happens to be equivalent to £7 [CISD, full]— £7[SCF] for Li, but is equal to £7 [CISD, full] - 
E [CISD, valence] for Be. 

There exists a tradition of optimizing basis sets for correlated calculations at the CISD 
level, since CISD is a variational method. However, optimizations for the CCSD or CCSD(T) 
counterparts of the abovementioned correlation energies yield very similar basis sets, which 
would have been of identical quality in molecular calculations. (A table of these exponents 
for the CVDZ, CVTZ, and CVQZ cases can be found in the Supplementary Material 39^.) 

By analogy with the cc- pCVnZ basis sets for B-NeQ, 1 sip, 2s2pld, 3s3p2dlf, and 
AsAp3d2flg sets of primitives were added to the cc-pVDZ, cc-pVTZ, cc-pVQZ, and cc-pV5Z 
basis sets, respectively. Successive angular momenta were optimized individually, after which 
the complete set of exponents was further optimized together. 

Multiple local minima exist for the larger sets, and care was taken to ensure that the 
final optimized basis set reflects the most 'contiguous' (or most 'even-tempered') solution in 
the sense of having no obvious 'gaps' between exponents. 

Aside from this issue, the basis set optimizations proceeded uneventfully. 



C. Na and Mg 

Likewise, for Na and Mg, unpublished cc-pVnZ basis sets of Feller j^sj were taken as 
the starting point. In this case, we optimized for the 2s2p correlation energy, found as 
£[CISD,2s2p3s]-£[CISD,3s] for Mg and as simply £[CISD,2s2p3s]-£[SCF] for Na. Note that 
the Is orbitals are constrained to be doubly occupied throughout: not only will the effect 
of this 'deep core' orbital on chemically significant properties be negligible, but including 
it in the optimization would bias the exponents towards the large — but chemically quite 
invariant — core-core correlation energy rather than the chemically more significant 2s2p 
core-core and core-valence correlation energies. 

Since the highest core correlation orbital is of p symmetry in this case, the basis sets 
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are rather larger, involving addition of lslpld, 2s2p2dlf, 3s3p3d2flg, 4s4p4:d3f2glh core 
correlation sets to the cc-pVDZ, cc-pVTZ, cc-pVQZ, and cc-pV5Z sets, respectively. 

For the CV5Z basis sets for Na and Mg, simply optimizing a set of core correlation s 
functions leads to such serious near-linear dependence problems that the basis would in 
practice be unusable. For want of an alternative, we simply uncontracted an additional four 
s primitives instead. 



D. K and Ca 



The core- valence gap in Ca, and particularly K, is so small that the optimization of regular 
cc-pVnZ basis sets would be of largely academic interest. Feller and coworkers previously 
published CVDZ, CVTZ, and CVQZ basis sets for kJ^; the sp exponents for the CVDZ and 
CVTZ basis sets were taken from (15s, 12p) and (18s, 15p) basis sets, respectively, of Ahlrichs 



and coworkers |42j], while those for the CVQZ basis set were taken from the Partridge-1 
set 3|, which is of (23s, 19p) uncontracted size. 

Initially we employed these basis sets unaltered, and merely added a CV5Z basis set. 
This was obtained as follows, starting from the Partridge-3 basis set: 

• the basis set was contracted to miminal (i.e. [4s3p]) in an atomic calculation on the 
ground state 

• the four outermost s and p primitives each were decontracted 

• [Ad3f2glh] functions were optimized at the valence-only CISD level for the K2 molecule 
at its experimental bond distance, yielding an intermediate cc-pV5Z basis set 

• core correlation functions were obtained by optimizing the atomic i?[CISD,3s3p4s]- 
£/[SCF] energy. The optimization of the s and p functions was plagued by insur- 
mountable near-linear dependence problems, and we ended up merely decontracting 
an additional four primitives of s and p symmetries each. In addition, a crossover 
occurred between the highest d exponent of the underlying cc-pV5Z basis set and the 
lowest d exponent required for the inner-shell correlation, and it was decided to merely 
keep the outermost three d primitives constant while optimizing all five remaining d 
primitives for 3s3p correlation. 
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Thus, our original CV5Z basis set was obtained. However, as evidenced in Table HJ 
basis set convergence of properties for K 2 and KH was less than satisfactory. In particular, 
inspection of the dissociation energy as a function of the basis set reveals that for K2, the 
CVTZ, CVQZ, and CV5Z results are nearly collinear, which is obviously undesirable. For 
this reason, the d, f, . . . exponents of the K CVTZ and CVQZ basis sets were reoptimized 
in the same manner as the CV5Z basis set. In the process we found the 'valence correlation' 
d and / exponents in the original Feller basis sets to be excessively 'tight', and with the 
adjusted basis sets (denoted 'version 0.1' in Table |l}, increases of 6.5 and 7.0 kJ/mol are 
seen in the CVTZ and CVQZ Dq values. Other molecular constants are also affected quite 
significantly. While the Feller CVTZ and CVQZ were not found to be as problematic for 
the other K compounds considered in this work, they are clearly unsuitable for K2 or any 
other system with similar long-distance bonding. 

Analogously, we optimized CVnZ (n=D,T,Q,5) basis sets for Ca, starting from Ahlrichs 
14s9p, Partridge-1 20sl2p, Partridge-2 23sl6p, and Partridge-3 26sl6p primitive sets. Ini- 
tial attempts to optimize basis sets for the 3 P atomic excited state (analogously to the 
sp-hybridized state of Be) yielded exceedingly poor basis set convergence of molecular prop- 
erties, and these basis sets were abandoned in favor of 1 S ground-state optimized basis sets. 
(We also attempted optimizations in the CaH 2 molecule, and found basis set parameters 
obtained there to be very similar to those for the 1 S atomic ground state.) Valence correla- 
tion functions were optimized for the valence correlation energy of the ground-state calcium 
atom, and core- valence functions for £'[CISD,3s3p4s]-£ , [CISD,4s]. For the CVQZ and CV5Z 
basis sets, the optimized 'valence' and 'subvalence' d and / shells interlock, and as a result 
only the outermost 2d2f were held at their valence-optimized values and the remainder 
reoptimized for inner-shell correlation. No such issues arose with the g and h functions. 
Once again, basis set convergence of molecular properties (e.g. for CaH, but also for other 
diatomics, not displayed here) was found to be somewhat unsatisfactory. 



For both atoms, we then optimized a sequence of four 'stretch-tuned' basis setsj4jj, 
i.e. basis sets in which the exponents of each angular momentum obey the following four- 
parameter relation: 

lnC„ = a + 7i(/3 + (n-l)( 7 + (n-2)<5)) (1) 

(A total of only eight parameters thus had to be optimized for each basis set.) The particular 
sequence of contraction sizes chosen was (20sl2p), (22sl4p), (24sl6p), and (2Qsl8p). The 



valence and inner-shell correlation functions were optimized as before. 

These basis sets (denoted 'version 0.2' in Table U} appeared to have better convergence 
properties. However, something was clearly still not satisfactory. For instance, basis set 
convergence for the ionization potentials of K and Ca atoms was found to be atypically slow 
(Table HH). 

It then occurred to us that, as the 4s orbital of these atoms is very diffuse, the atom- 
optimized p functions may be too 'tight' (as they were optimized for atoms with vacant 
4p orbitals) to contain suitable first polarization/angular correlation functions for the outer 
part of the 4s orbital. This assumption was verified by optimizing single 'probe' basis 
functions of various symmetries on top of CVnZ basis sets held constant. (Only valence 
electrons were correlated, as otherwise the functions would have fallen into the 'gravity 
well' of the inner-shell electrons.) A very noticeable lowering of the total energy was seen 
upon adding a single p function, the optimum exponent of which indeed turned out to be 
'looser' than any of the existing p functions. An additional 'loose' p function led to a further 
improvement for the K CVnZ (n=D,T,Q) basis sets, but not for the K CV5Z basis set (the 
outer p primitive of which is already quite 'loose') or any of the Ca basis sets. Consequently, 
the valence correlation parts of the CVnZ basis sets were reoptimized with two and one 
additional p functions added, respectively, for K and Ca, and afterwards the inner-shell part 
was reoptimized for consistency. Thus our final CVnZ basis sets were obtained. As seen in 
Table HI basis set convergence of molecular parameters is now satisfactory. For IP(K) and 
IP(Ca), nearly exact values are now obtained (Table ITT]) . 

Since the primitive K and Ca basis sets may be somewhat bulky for some applications, 

we have also generated SDB-CVnZ basis sets, which are combinations of the valence and 

I — I 

subvalence parts of the CVnZ basis sets with small-core Stuttgart-Dresden-Bonn |45| rel- 
ativistic effective core potentials. The latter replace the (Is2s2p) electrons. Unlike the 
situation with large-core basis set where reoptimization of the basis set is essential 22], we 
have kept all exponents at their all-electron values and merely 'pruned' and recontracted the 
basis set. Specifically, the atomic SCF calculation was repeated with the (Is2s2p) orbitals 
replaced by the SDB pseudopotential denoted 'ECP10MWB'. Then all primitives with ab- 
solute coefficients of less than 10~ 5 in any valence orbitals were deleted, and the orbitals 
recontracted. 

All basis sets obtained in this work are available in machine-readable form (Gaussian and 
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MOLPRO formats) as supplementary material |39| to the present paper. Contracted basis 
set sizes and numbers of basis functions for each element are given in Table IIV1 



IV. VALIDATION FOR DIATOMIC MOLECULES 
A. Diatomic metal hydrides 

Since hydrogen has no core electrons, there are no differential core-core contributions 
(just core-valence) to the molecular properties. Computed and experimental data for the 
diatomic hydrides are given in Table IVl 

As noted previously, the inner-shell contribution to the LiH spectroscopic constants is 
quite appreciable, as can be expected from the small Li(ls)-H(2s) gap of about 2 a.u. We 
note that the CCSD(T,riv)/CV5Z results reproduce the Born-Oppenheimer bond distance 
and harmonic frequency to within 0.0002 A and 0.2 cm -1 , respectively. Neglect of Li(ls) cor- 
relation leads to errors of 0.012 A and 13.5 cm -1 , respectively. Nevertheless, the contribution 
to D e does not exceed 1.2 kJ/mol. 

The corresponding core-valence gap for BeH is almost doubled, which translates into sub- 
stantially reduced core- valence contributions (relatively speaking) to r e and u e . Interestingly, 
the contribution to Do reaches nearly 2 kJ/mol. 

In contrast, the Na(2p)-H(ls) gap narrows to no more than 1 a.u., and in NaH we see 
inner-shell correlation contributions of 33 cm -1 to u> e and 0.035 A to r e . Interestingly, once 
more D e is nearly unaffected. Once more, the CCSD(T,riv)/CV5Z results are in excellent 
agreement with experiment for the spectroscopic constants: applying a W2-type extrapola- 
tion for Dq results in a value of 182.35 kJ/mol, within 0.3 kJ/mol of experiment. 

Given that the Mg(2p)-H(ls) gap is not dissimilar from the corresponding one in LiH, 
it is not greatly surprising that the importance of core-valence contributions is again much 
smaller than in NaH. However, the 8 kJ/mol contribution to D is considerably more sig- 
nificant. 

In KH, the K(3p) and H(ls) orbitals are closer than 0.5 a.u., and the contributions of 
nearly 0.1 A and 44 cm -1 to r e and u) e , respectively, clearly suggest that any sort of 'valence 
correlation only' calculation on a K compound should be viewed with great suspicion. Note 
that the core-valence contribution to Dq is negative in this case, as was previously found 
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(e.g. 12]) for aluminum and silicon hydrides. 

In comparison, the Ca(3p)-H(ls) gap widens to 0.84 a.u., and core- valence effects on the 
molecular properties are mitigated accordingly. Like for MgH and KH, inner-shell correlation 
in fact slightly reduces the dissociation energy. 

With the exception of KH and CaH, the differential core-valence contributions appear 
to have converged with respect to the basis set at the CVQZ stage. This in itself satis- 
fies a minimum requirement for treating the inner-shell correlation contribution separately, 
something which is definitely inappropriate for KH. 

On the whole, agreement between CCSD(T,riv)/CV5Z and experiment can only be de- 
scribed as quite satisfactory. For systems other than LiH and NaH, imperfections in the 
CCSD(T) electron correlation method may account for most of the remaining discrepancy 



between computed and observed harmonic frequencies 



46|. 



B. Metal diatomics 



The alkaline earth metal dimers (particularly Be2) exhibit such severe multireference 
character that they warrant studies in themselves (which would, however, focus on electron 
correlation methods rather than basis sets, see e.g. j^J and references therein). Results for 
species other than Be2 can be found in Table IVI1 

In LI2, subvalence correlation accounts for 0.024 A, which is sizable by spectroscopist's 
standards. In Na2, this is drastically increased to nearly 0.1 A, and reaches a whopping 
0.23 A in K 2 . Interestingly, these contributions are quite close to double their metal hydride 
counterparts. As expected, LiNa and NaK represent scenarios intermediate between the 
homonuclear diatomics of the constituent elements. Changes in u e are modest in absolute 
numbers, but for these molecules represent relative errors of up to 5%. It is quite clear how- 
ever that, with the possible exception of Li2, inclusion of subvalence correlation is essential 
for reliable molecular parameters. Once again we see, however, that bond energies are only 
affected mildly: 0.8 kJ/mol in Li 2 , 1.7 kJ/mol in LiNa, 1.6 kJ/mol in Na2, and only 0.5 
kJ/mol in K 2 . 

Once again, a basis set of CVQZ quality appears to be close to convergence for the 
differential core-valence effects. 

Effects for Mg2 and Ca2 are considerably milder than for their alkali neighbors, yet (3s, 3p) 
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correlation reduces the Ca 2 bond distance by 0.07-0.08 A, not negligible by any reasonable 
standard. The discrepancies between theory and experiment for Mg 2 primarily reflect the 
inadequacy of the CCSD(T) method for this highly multireference system. In contrast, Ca 2 
yields quite satisfying results compared to experiment. 

Agreement between our best calculations and experiment can only be described as excel- 
lent for Li 2 . The same would be true of Na 2 if it were not for the bond length which is still 
0.003 A too long at the CCSD(T,riv)/CV5Z level. However, at the DK-CCSD(T)/CVQZ 
level, we find a scalar relativistic correction to the bond length of -0.006 A. (If this correc- 
tion would seem to be exaggerated, we note that the Na 2 potential curve is quite flat, and 
that already for OH - and HF, scalar relativistic corrections were found to be required for 
spectroscopic-quality results 43].) In K 2 , our best calculated results are likewise in excellent 
agreement with the experimental spectroscopic data except for r e , which is calculated to be 
0.014 A longer than the observed value. At the DK-CCSD(T)/CVQZ level, we find a rela- 
tivistic contribution of -0.013 A to the bond length, which explains most of the discrepancy. 
In addition, the (2s, 2p) 'deep core' orbitals are energetically in the same range as the (Is) 
orbitals in C and N, so it cannot be entirely ruled out that 'deep core' correlation may affect 
molecular properties in K compounds. 

We optimized a CCVTZ (deep core valence triple zeta) basis set for K. In these calcu- 
lations, the CVTZ basis set was held constant and exponents of 2s2p2dlf basis functions 
optimized for £'[CISD,2s2p3s3p4s] — £7[CISD,3s3p4s], i.e. for deep-core correlation energy 
taken in isolation. As can be seen in Table ITPTJ these exponents are obviously much 'tighter', 
by almost an order of magnitude, than those required for subvalence correlation. Hence, any 
'deep core' result in a mere valence and subvalence CVnZ (let alone a valence-only cc-pVnZ 
set) should be regarded with skepticism at best. 

The only change of note we see in the molecular properties when using this CCVTZ 
basis set for K 2 is a shortening of the K K bond by another 0.0018 A. We also attempted 
to optimize a CCVQZ basis set but ran into insurmountable near-linear- dependence prob- 
lems; however, the CCVTZ result should at least give an indication of the (fairly modest) 
magnitude of this effect. 
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C. Diatomic metal halides 

This set of molecules is particularly revelant because of the deep-lying valence s orbitals 
on the halogen atoms. Relevant results are collected in Table IVHI 

Firstly, for LiF the contribution of Li (Is) correlation to the binding energy is somewhat 
more noticeable (4 kJ/mol). Its inclusion also causes a contraction of the bond by 0.015 A 
and an increase in u e of 12 cm -1 . The further effect of permitting F(ls) correlation is nearly 
negligible in comparison. The Cl(2s,2p) correlation effect is somewhat more noticeable but 
still Li(ls) correlation accounts for the lion's share of the changes. 

The F(2s) and Na(2p) orbitals are nearly degenerate, and hence one would expect valence- 
only calculations on NaF to be problematic, to say the least. Nevertheless, while the effect of 
Na(2s,2p) correlation on the molecular constants is quite noticeable and its inclusion clearly 
essential for accurate calculations, the results obtained are not outright nonsensical. Of 
course, considering the difference in electronegativity between the elements, the combination 
of Na + and F~ would yield a more appropriate (if perhaps extreme) 'atoms in molecules' 
picture — and the corresponding orbital energy gap in this admittedly extreme scenario is 
0.72 a.u. with the Partridge 3 basis set. The corresponding 'ionic' gap for NaCl amounts to 
1.07 a.u., explaining why also for NaCl, valence-only results are not totally unreasonable. 

The person carrying out calculations on KF and KC1 has no such luck. The inclusion 
of K(3s3p) correlation shifts u e in KF upwards by no less than 20%, and the contribution 
to D e is seen to be on the order of 170 kJ/mol. In fact, the 'valence' results are better 
described as referring to correlating the z component of K(3p) and freezing the F(2s), which 
obviously makes no chemical sense at all. The valence-computed dissociation energy of KC1 
agrees deceptively well with experiment, but an error of 0.3 A in the bond distance and of a 
factor of three in the anharmonicity constant should discourage any quantum chemist from 
performing valence-only calculations on this type of species. 

While obviously significant for spectroscopic purposes, the effects of metal subvalence 
correlation in MgF and MgCl are clearly less prominent than in NaF and NaCl, respectively, 
and with the beryllium halides, one clearly could 'get away' with a differential treatment. 
In contrast, in CaF a calcium subvalence contribution to D of 180 kJ/mol is found, not to 
mention +46 cm -1 on u e and a factor of two on u e x e . Oddly, the effect on the bond length 
is comparatively small (0.04 A). The opposite scenario is seen for CaCl, with a Ar e of 0.1 A 
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but otherwise not outlandishly erroneous values for the other spectroscopic constants. CaF 
and CaCl results resemble more the potassium halides. 

And once again, the CCSD(T,riv)/CV5Z results agree as well with experiment as can 
reasonably be expected. 

We have also considered calculations in which subvalence correlation on the halide was 
permitted. However, for K and Ca compounds with {C, N, O, F}, the (Is) orbital of these 
latter elements is in fact below the K(2s,2p) 'deep core' orbitals in energy, and therefore any 
such calculation will require the use of CCVnZ 'deep core' basis sets on the metal. Such 
calculations are sufficiently costly that one might want to limit them to single-point energy 
calculations in high-accuracy thermodynamic work. For MCI, the most noticeable effect of 
including (2s, 2p) subvalence correlation on chlorine is a contraction of r e by 0.003 A. 



D. The metal chalcogenides 

We will primarily focus on the oxides. Results can be found in Table IVIIII 
Clearly, changes of 13 cm -1 and 6 kJ/mol, respectively, for u e and D e of BeO suggest 
the importance of metal (Is) correlation in this system; contributions from O(ls) correla- 
tion are an order of magnitude less important. Note that the experimental data can be 
reproduced almost exactly at the CCSD(T,all)/CV5Z level, despite BeO having pronounced 
multireference character and in fact being on the borderline of applicability of CCSD(T) (see 
e.g. 4||). In the isovalent BeS system, Be(ls) correlation is still noticeably more important 
than S(2s,2p) correlation, although the difference is not as pronounced as its counterpart in 
BeO. 

Changes between valence-only and RIV treatment of 0.017 A in r e (MgO) and 23 cm -1 
in (jj e (MgO) speak for themselves. The MgO molecule has very pronounced multireference 
character, and accurate reproduction of the experimental spectroscopic constants would 
required an elaborate multireference treatment. Subvalence correlation effects in MgS are 
rather milder. 

In contrast, inspection of the computed spectroscopic constants for CaO immediately 
reveals that valence-only calculations on this system are essentially exercises in wasting 
computer time. Valence-only bond distances and dissociation energies are off by a quarter 
of an A, and a factor of six in Dq. Respectable agreement with experiment can however 
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be reached by correlating the Ca(3s,3p) electrons. Obviously, the large uncertainties in the 
measured dissociation energies preclude really fine comparisons with experiment: on the 
basis of the preceding, it can perhaps be stated that the calculated binding energies are 
more reliable than the experimental ones for the alkali(ne earth) metal chalcogenides. 

Once again, correlating the (Is) orbital on O would require including (2s, 2p) deep-core 
correlation in Ca, otherwise meaningless results are obtained. 



E. Additional d functions on Ca 



Wesolowski et al.[50], in a study of calcium oxide, observed especially poor performance 



for a valence-only atomic natural orbital (ANO) basis set This was considerably im- 
proved when high-exponent d functions were added. They noted that Ca + has a low-lying 
2 D state only 1.7 eV above the 2 S ground state, and showed that said ANO basis set over- 
estimates this separation by nearly a factor of three. Addition of several very high-exponent 
d functions, taken from the 6-311G(2df) basis set for Ca 52] (which treats the (3d) orbital 
on a valence footing), results in an at least semiquantitatively correct result. 

As the present basis sets were developed from the ground up with inner-shell correlation 
in mind, they intrinsically contain high-exponent d functions and should therefore be much 
less susceptible to the problem. As shown in Table IIX[ quite good agreement with the 
experimental transition energy is obtained for CVQZ and CV5Z basis sets, but the smaller 
basis sets leave something to be desired. 

We therefore proceeded to optimize CV(n + d)Z, CV(n + 2<f)Z, and CV(n + 3d)Z basis sets 
in the following fashion. (Both the nomenclature and the procedure bear resemblance to 

n 

the cc-pV(n + d)Z basis sets for second-row elements [22J.) All angular momenta other than 
d were kept constant, as were the valence-optimized d functions (see above). The remaining 
d functions were reoptimized for inner-shell correlation with one, two, or three additional d 
functions added. (In the CV(Q+d)Z and CV(Q+2d)Z basis sets, it was found necessary to 
constrain the exponents to stretch-tuned and even-tempered sequences, respectively.) 

As can be seen in Table El the 2 D *- 2 S excitation energy indeed converges considerably 
faster for the CV(n + d)Z and CV(n + 2d)Z series than for the CVnZ series. After applying 
a scalar relativistic correction by means of the Douglas-Kroll method in an uncontracted 
CV(Q+d)Z basis sets, very good agreement with experiment can be achieved. As for the 
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molecular properties of CaO fTable |VHI|) . the effect of additional d functions is quite dra- 
matic in the CVDZ case and still quite important in the CVTZ case, but tapers off for CVQZ 
and becomes quite insignificant for CV5Z. We noticed similar behavior in the other CaX 
systems considered in this paper (Tables W \ Wl \ and IVIIj) . Convergence as a function of the 
number of additional d functions appears to be approached for CV(D+3d)Z, CV(T+2d)Z, 
CV(Q+d)Z, and CV5Z (note particularly Table IVrTT]) . Basis set convergence of molecular 
properties along this sequence is clearly much smoother than for the unmodified CVnZ basis 
sets. This is particularly true for CaO and, to a somewhat lesser extent, for CaF, least so 
for CaH and Ca2 where basis set convergence was adequate to begin with. 

V. CONCLUSIONS 

Core-valence basis sets have been developed for the alkali and alkaline earth metals Li, 
Be, Na, Mg, K, and Ca. Validation calculations for a number of diatomics involving these 
systems reveal basis set convergence to be satisfactory. 

The addition of low-exponent p functions to the K and Ca basis set is found to be essential 
for smooth basis set convergence of molecular properties. These functions accommodate 
angular correlation from the outer part of the 4s orbital. 

Valence-only calculations on K and Ca chalcogenides and halides yield large errors at best, 
and chemically nonsensical results at worst. In general, inclusion of subvalence correlation 
in K and Ca compounds is absolutely essential for even reliable (let alone accurate) results, 
while for accurate calculations, subvalence correlation in Na (and, to a lesser extent, Li 
and Mg) should be included, preferably as a 'baseline' treatment rather than an additive 
correction. The latter appears to be sufficient for Be compounds. Any study that seeks to 
address the importance of first-row element core correlation in K and Ca compounds should 
also consider deep-core (2s,2p) correlation from these latter elements. 

For compounds involving Ca bound to highly electronegative elements, we strongly rec- 
ommend basis sets with additional d functions on Ca: our suggested sequence for basis set 
convergence studies would be CV(D+3d)Z, CV(T+2d)Z, CV(Q+d)Z, CV5Z. 
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TABLE I: Comparison of different CVnZ basis sets for properties of K2 and KH. 







r e 




LU e X e 


D 






(A) 


(cm ) 


(cm" 1 ) 


(kj/mol) 


KH 


CVDZ Feller 




2.2835 


935.7 


13.80 


124.3 


CVDZ version 0.1 




2.2847 


930.0 


13.87 


123.2 


CVDZ version 0.2 




2.2673 


960.5 


18.61 


125.4 


Final CVDZ (incl. 


+2p) 


2.3232 


902.1 


11.42 


136.1 


CVTZ Feller 




2.2382 


998.8 


16.69 


156.7 


CVTZ version 0.1 




2.2410 


991.0 


17.11 


153.0 


CVTZ version 0.2 




2.2360 


1025.3 


13.64 


159.5 


Final CVTZ (incl. 


+2p) 


2.2591 


976.3 


15.41 


161.2 


CVQZ Feller 




2.2353 


1002.9 


14.24 


166.0 


CVQZ version 0.1 




2.2353 


1003.5 


14.73 


165.0 


CVQZ version 0.2 




2.2419 


991.7 


13.93 


166.2 


Final CVQZ (incl. 


+2p) 


2.2479 


979.0 


14.32 


167.8 


CV5Z version 0.1 




2.2406 


989.8 


15.31 


168.8 


CV5Z version 0.2 




2.2445 


985.3 


15.81 


169.4 


Final CV5Z (incl. 


+2p) 


2.2461 


982.6 


15.84 


170.0 


Expt. a 




2.240164(10) 986.6484(41) 


15.54 


170.972 6 


K 2 


CVDZ Feller 




4.1250 


77.1 


0.45 


27.5 


CVDZ version 0.1 




4.1507 


74.6 


0.42 


27.2 


CVDZ version 0.2 




4.2223 


69.9 


0.43 


24.5 


Final CVDZ (incl. 


+2p) 


4.0511 


88.5 


0.29 


50.0 


CVTZ Feller 




3.9956 


87.2 


0.40 


36.9 


CVTZ version 0.1 




3.9561 


92.1 


0.32 


43.4 


CVTZ version 0.2 




4.0180 


78.1 


0.38 


30.9 


Final CVTZ (incl. 


+2p) 


3.9837 


90.1 


0.25 


50.2 


CVQZ Feller 




3.9434 


91.8 


0.34 


41.7 


CVQZ version 0.1 




3.9522 


93.0 


0.29 


48.7 


CVQZ version 0.2 




3.9502 


90.9 


0.35 


44.0 


Final CVQZ (incl. 


+2p) 


3.9557 


92.0 


0.33 


51.2 


CV5Z version 0.1 




3.9440 


92.9 


0.29 


51.6 


Final CV5Z (incl. 


+2p) 


3.9380 


93.0x 


0.34 


52.2 



Expt. c 3.92435(2) 92.39766(47) 0.32485(12) 53.243(2) 

(a) Uehara, H., Horiai, K., and Konno, T., 1997, J. Mol. Struct. 413, 457. 

(b) Ref . l55l . corrected to K using spectroscopic constants given to the left. 

(c) Amiot, C, Verges, J., and Fellows, C. E., 1995, J. Chem. Phys. 103, 3350. 
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TABLE II: Basis set convergence of ionization potentials of K and Ca (eV), with and without 
low-exponent p functions added. 



K: IP(expt.)=4.34077 a 





CCSD(T) 


DK-CCSD(T) 






+ 2p 




+2p 


CVDZ 


4.068 


4.168 


4.086 


4.185 


CVTZ 


4.218 


4.275 


4.236 


4.292 


CVQZ 


4.300 


4.301 


4.317 


4.318 


CV5Z 


4.308 


4.308 


4.326 


4.326 


CVooZ 


4.316 


4.315 


4.335 


4.334 



Ca: IP(expt.)=6.11316 a 





CCSD(T) 


DK-CCSD(T) 






+ IP 




+lp 


CVDZ 


5.598 


6.003 


5.619 


6.022 


CVTZ 


5.797 


6.050 


5.818 


6.069 


CVQZ 


6.038 


6.075 


6.085 


6.094 


CV5Z 


6.073 


6.083 


6.094 


6.102 


CVooZ 


6.110 


6.091 


6.103 


6.110 



Basis set limits using expression |56l : E x El + {El — El-i)/{{L/L — l) 3 — 1) 
(a) Handbook of Chemistry and Physics (CRC Press, Boca Raton, FL, 2001), p. 10-175. 



TABLE III: Valence correlation, core-valence and deep-core correlation exponents in the CVTZ 
and CCVTZ basis sets for potassium. 





Valence 


Core 


Deep core 


s 


0.0181, 0.0501 


0.346, 1.779 


4.281, 20.055 


p 


0.0275, 0.0598 


0.529, 2.599 


7.240, 16.761 


d 


0.0599, 0.431 


1.036, 2.586 


9.724, 28.720 


f 


0.0888 


1.129 


15.150 



TABLE IV: Contracted basis set sizes and total numbers of basis functions for the CVnZ basis 
sets for group 1 and 2 elements 





CVDZ 


CVTZ 


CVQZ 




CV5Z 




Li,Be 


[4s3pld] 18 


[6s5p3dlf] 43 


[8s7p5d3flg] 


84 


[10s9p7d5f3glh] 


145 


Na,Mg 


[5s4p2d] 27 


[7s6p4d2f] 59 


[9s8p6d4f2g] 


109 


[Ilsl0p8d6f4g2h] 


181 


K,Ca 


[6s5p2d] 31 


[8s7p4d2f] 63 


[10s9p6d4f2g] 


113 


[12sllp8d6f4g2h] 


185 
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TABLE V: Computed and observed spectroscopic constants (kJ/mol, A, cm 1 as appropriate) for 
alkali and alkaline earth metal hydrides. 

















Do 




Val RIV 


Val 


RIV 


Val 


RIV 


Val 


RIV 


LiH 


CVDZ 


1.6191 1.6106 


1369.1 


1379.6 


20.80 


21.55 


209.1 


210.6 


CVTZ 


1.6081 1.5994 


1394.9 


1400.1 


22.72 


22.08 


227.6 


228.7 


CVQZ 


1.6069 1.5960 


1393.1 


1405.3 


23.61 


22.84 


231.4 


232.6 


CV5Z 


1.6074 1.5954 


1391.8 


1405.7 


22.78 


23.21 


232.4 


233.6 


Expt. a 


1.5955991(1) 


1405.49805(76) 


23.167899(714) 


244.33 


, 234.354(4)9 


BeH 


CVDZ 


1.3587 1.3556 


2036.1 


2033.9 


34.35 


34.26 


177.4 


178.1 


CVTZ 


1.3498 1.3448 


2039.9 


2053.3 


35.40 


35.79 


191.3 


193.1 


CVQZ 


1.3458 1.3411 


2054.9 


2065.6 


36.32 


36.65 


195.9 


197.7 


CV5Z 


1.3456 1.3407 


2055.3 


2065.9 


36.57 


36.68 


196.8 


198.7 


Expt. 6 


1.3411(1) 


2061.66 


37.18 


196.2, 


193.3 ± 1.39 


NaH 


CVDZ 


1.9237 1.9078 


1116.3 


1115.7 


17.41 


17.58 


158.0 


158.7 


CVTZ 


1.9232 1.8949 


1131.6 


1154.0 


16.46 


18.56 


176.1 


176.5 


CVQZ 


1.9256 1.8907 


1122.8 


1159.3 


10.10 


17.89 


179.7 


180.2 


CV5Z 


1.9235 1.8883 


1137.8 


1171.1 


18.43 


19.67 


181.7 


181.2 


Expt. c 


1.88703 


1171.968 ±0.012 


19.703 ±0.010 


(181), 182.03 ±0.25 9 


MgH 


CVDZ 


1.7442 1.7398 


1468.1 


1458.9 


27.22 


26.31 


107.6 


106.3 


CVTZ 


1.7434 1.7332 


1496.6 


1499.7 


28.36 


28.59 


123.3 


120.3 


CVQZ 


1.7413 1.7299 


1494.5 


1498.1 


28.46 


29.10 


127.1 


122.9 


CV5Z 


1.7401 1.7289 


1493.2 


1498.8 


28.25 


29.31 


128.3 


123.5 


Expt. d 


1.729828(2) 


1495.2632(22) 


31.64139(84) 


129, 


122.7 ± 2. 9 9 


KH 


CVDZ 


2.3567 2.3232 


907.6 


902.1 


12.58 


11.42 


141.0 


136.1 


CVTZ 


2.3230 2.2591 


950.7 


976.3 


14.97 


15.41 


167.5 


161.2 


CVQZ 


2.3241 2.2479 


944.4 


979.0 


14.41 


14.32 


173.3 


167.8 


DK-CVQZ 


2.2447 




981.7 




14.38 




167.4 


CV5Z 


2.3274 2.2461 


938.8 


982.6 


14.81 


15.84 


174.8 


170.0 


Expt. e 


2.240164(10) 


986.6484(41) 


15.54 


179, 170.972s 



(continued on next page) 
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TABLE V: (continued) 



CaH 



CVDZ 


2.1453 


2.0941 


1231.5 1227.4 


20.04 


18.17 


114.3 


115.7 


CV(D+d)Z 




2.0539 


1249.6 




18.21 




121.8 


CV(D+2d)Z 




2.0323 


1260.3 




18.25 




126.7 


CV(D+3d)Z 




2.0227 


1290.5 




20.56 




129.8 


CVTZ 


2.1341 


2.0415 


1253.2 1273.9 


18.70 


17.29 


141.5 


148.4 


CV(T+d)Z 




2.0204 


1279.1 




16.85 




153.7 


CV(T+2d)Z 




2.0149 


1283.8 




17.46 




154.7 


CVQZ 


2.0899 


2.0063 


1270.4 1297.9 


20.62 


18.80 


158.5 


162.5 


CV(Q+d)Z 




2.0048 


1298.7 




18.72 




162.9 


DK-CVQZ 




2.0077 


1296.9 




18.53 




160.4 


DK-CV(Q+d)Z 




2.0062 


1297.9 




18.65 




160.8 


CV5Z 


2.0536 


2.0027 


1271.9 1299.1 


18.68 


18.93 


169.3 


165.2 


ExptJ 


2.002366(16) 


1298.3999(40) 19.1842(28) 


< 164, 164.19 



(a) (LiH) Dulick, M., Zhang, K.-Q., Guo, B., and Bernath, P.F., 1998, J. Mol. Spectrosc. 188, 14. 

(b) (BoH) (1) Colin, R. and De Greef, D., 1975, Can. J. Phys. 53, 2142. (2) Martin, J. M. L., 1998, Chem. Phys. Lett. 283, 
283. 

(c) (NaH) Pesl, F. P., Lutz, S., and Bergmann, K., 2000, Ever. Phys. J. D. 10, 247. 

(d) (MgH) Lemoine, B., Demuynck, C, Destombes, J.L., and Davies, P.B., 1988, J. Chem. Phys. 89 , 673. 

(e) (KH) Uehara, H., Horiai, K., and Konno, T., 1997, J. Mol. Struct. 413, 457. 

(f) (CaH) Petitprez, D., Lemoine, B., De muy nck, C, Destombes, J.L., and Macke, B., 1989, J. Chem. Phys. 91, 4462. 

(g) First value from Huber and Herzberg^J, second from CRC Handbook lssl corrected to K with expt. ui e a> e and oj e . 



TABLE VI: Computed and observed spectroscopic constants (kJ/mol, A, cm 1 as appropriate) for 
alkali and alkaline earth metal diatomics. 









e 






Do 




Val RIV 


Val 


RIV 


Val 


RIV 


Val RIV 




CVDZ 


2.7284 2.6963 


340.4 


346.7 


2.32 


2.44 


92.0 95.5 


CVTZ 


2.7006 2.6802 


344.7 


348.8 


2.50 


2.53 


97.3 98.2 


CVQZ 


2.6984 2.6756 


346.6 


351.0 


2.53 


2.58 


98.5 99.2 


CV5Z 


2.6985 2.6741 


346.6 


351.3 


2.54 


2.59 


98.8 99.4 


Expt. a 


2.6729 


351.43 


2.610 


100.9, 107.45 ±A h 


LiNa 


CVDZ 


2.9679 2.9386 


244.1 


247.5 


1.38 


1.45 


77.0 79.9 


CVTZ 


2.9505 2.9051 


247.6 


253.5 


1.50 


1.55 


81.0 82.5 


CVQZ 


2.9496 2.8953 


248.6 


256.1 


1.56 


1.64 


81.7 83.0 


CV5Z 


2.9500 2.8910 


248.9 


256.4 


1.55 


1.61 


82.0 83.1 


Expt. 6 


2.88881(2) 


256.4577(14) 


1.5808(6) 


N/A, 84.728 ± 0.001 h 



Na 2 



CVDZ 


3.2048 3.1829 


149.9 


152.8 


0.80 


1.48 


65.0 67.9 


CVTZ 


3.1779 3.1007 


151.7 


158.4 


0.66 


0.69 


68.4 70.9 


CVQZ 


3.1773 3.0884 


151.7 


158.7 


0.69 


0.74 


68.8 70.6 


DK-CVQZ 


3.0819 




159.4 




0.77 


70.9 


CV5Z 


3.1783 3.0822 


152.1 


159.3 


0.67 


0.71 


69.1 70.7 


Expt. c 


3.0795(1) 


159.103(3) 


0.7190(6) 


69, 71.0173 ±0.0001 h 


NaK 


CVDZ 


3.6793 3.6072 


115.7 


118.4 


0.45 


0.44 


58.0 59.6 


CVTZ 


3.6666 3.5379 


116.6 


121.3 


0.43 


0.44 


61.3 60.9 


CVQZ 


3.6671 3.5185 


116.7 


122.8 


0.44 


0.44 


61.8 61.4 


CV5Z 


3.6661 3.5070 


117.0 


124.0 


0.44 


0.44 


61.9 61.9 


Expt. d 


3.49903 


123.993 


0.3045 


60, 64.089 ± 0.008 h 








M 


g2 






CVDZ 


4.9450 4.8450 


17.1 


18.8 


0.60 


0.62 


1.0 1.3 


CVTZ 


4.0581 4.0333 


41.0 


43.6 


1.65 


1.51 


3.0 3.3 


CVQZ 


3.9748 3.9717 


46.6 


46.0 


1.68 


1.64 


3.9 3.8 


CV5Z 


3.9579 3.9709 


47.6 


45.9 


1.67 


1.50 


4.1 3.9 


Expt. e 


3.8905 


51.121 


1.645 


4.83, 7.037 ± 0.004' 1 


K 2 


CVDZ 


4.1626 4.0511 


85.4 


88.5 


0.28 


0.29 


48.7 50.0 


CVTZ/ 


4.1684 3.9837 


85.8 


90.1 


0.27 


0.25 


51.1 50.2 


CVQZ 


4.1694 3.9557 


85.8 


92.0 


0.28 


0.33 


51.6 51.2 


DK-CVQZ 


3.9407 




92.9 




0.21 


52.1 


CV5Z 


4.1693 3.9380 


85.6 


93.0 


0.28 


0.34 


51.7 52.2 


Expt.f 


3.92435(2) 


92.39766(47) 0.32485(12) 49.6, 52. 88(2)\ 53.243(2)9 



(continued on next page) 
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TABLE VI: (continued) 















D 






Val RJV 


Val 


RIV 


Val 


RIV 


Val 


RIV 


Ca 2 


CVDZ 


5.0096 4.7730 


26.8 


33.8 


1.11 


1.23 


2.5 


3.4 


CV(D+d)Z 


4.7956 




32.5 




1.17 




3.5 


CV(D+2d)Z 


4.7642 




33.5 




1.19 




3.7 


CV(D+3d)Z 


4.7554 




33.9 




1.17 




3.8 


CVTZ 


4.5409 4.4185 


53.6 


56.1 


1.05 


1.07 


8.3 


9.0 


CV(T+d)Z 


4.3932 




57.2 




1.08 




9.6 


CV(T+2d)Z 


4.3899 




57.2 




1.09 




9.6 


CVQZ 


4.4385 4.3440 


60.6 


60.9 


1.05 


1.07 


10.7 


11.0 


CV(Q+d)Z 


4.3417 




61.0 




1.07 




10.9 


DK-CVQZ 


4.3313 




59.1 




0.99 




11.0 


DK-CV(Q+d)Z 


4.3289 




59.2 




1.00 




11.1 


CV5Z 


4.4005 4.3261 


63.0 


61.5 


1.08 


1.07 


12.1 


11.5 


Expt.* 


4.2773, 4.277 


64.93, 65.07 


1.065, 1.09 


12.4±1, 13.1 



(a) (L12) Huber and Herzberg[^3] . 

(b) (LiNa) Fellows, C.E., 1991, J. Phys. Chem. 94, 5855. 

(c) (Na 2 ) Babaky, O. and Hussein, K. J., 1989, Can. J. Phys. 67, 912. 

(d) (NaK) Krou-Adohi, A. and Giraud-Cotton, S., 1998, J. Mol Spectrosc. 190, 171. 

(e) (Mg2) Huber and HerzberglsH, 

(f) (K2) With CCVTZ deep-core correlation basis set and correlating all but (Is) electrons (RIV values given in parentheses): 
r e =3.9809(3.9827) A, w e =90.2(90.2) cm" 1 , a; e x e =0.25(0.25) cm" 1 , D =50.3(50.2) kj/mol. (g) (K 2 ) Amiot, C, Verges, J., 
and Fellows, C.E., 1995, J. Chem. Phys. 103, 3350. 

(h) First value from Huber and Herzbergl57|. second from CRC Handbook 55J corrected to K with expt. u e x e and uj e . 

(i) First set of values: Huber and Herzberg 57], from Balfour, W. J. and Whitlock, R. F., 1975, Can. J. Phys. 53, 472; 
Second set of values: Vidal, C. R., 1980, J. Chem. Phys. 72, 1864, Bondybey, V. E. and English, J. H., 1984, Chem. Phys. 
Lett. Ill, 195. 
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TABLE VII: Computed and observed spectroscopic constants (kJ/mol, A, cm 1 as appropriate) 
for metal halides. 











UJ e X e 


D Q 






Val RIV All 


Val 


RIV All 


Val RIV All 


Val RIV 


All 


LiF 


CVDZ 


1.6074 1.5983 1.5977 


870.6 


876.1 876.9 


8.16 8.18 8.18 


539.7 543.8 


544.8 


CVTZ 


1.5903 1.5747 1.5741 


885.5 


891.2 891.5 


7.79 7.70 7.74 


558.7 562.4 


563.0 


CVQZ 


1.5817 1.5671 1.5665 


894.4 


906.1 906.6 


8.00 8.12 8.13 


568.9 572.7 


573.1 


CV5Z 


1.5808 1.5656 1.5650 


895.9 


908.0 908.6 


8.04 8.12 8.14 


571.0 575.0 


575.6 


Expt. a 


1.5638648(3) 


910.57272(10) 


8.207956(46) 


570, 537 ± 21 m 




LiCl 


CVDZ 


2.0900 2.0819 2.0806 


600.9 


607.4 608.7 


3.57 3.97 3.95 


434.5 436.3 


437.4 


CVTZ 


2.0486 2.0375 2.0351 


626.9 


633.1 634.3 


4.26 4.33 4.34 


454.3 455.9 


456.5 


CVQZ 


2.0430 2.0291 2.0266 


632.2 


640.1 641.0 


4.38 4.45 4.45 


465.5 467.4 


468.0 


Expt. ft 


2.0206719(2) 


642.95821(14) 


4.475085(57) 


467, 466 ± 13 m 




BeF 


CVDZ 


1.4148 1.4083 1.4078 


1190.9 


1197.6 1198.9 


8.27 8.52 8.54 


509.9 513.4 


514.6 


CVTZ 


1.3741 1.3696 1.3690 


1240.1 


1245.3 1246.4 


9.13 9.17 9.19 


546.2 549.2 


549.9 


CVQZ 


1.3689 1.3635 1.3531 


1254.8 


1264.1 1264.3 


9.30 9.39 9.40 


558.8 562.8 


563.1 


CV5Z 


1.3678 1.3623 1.3617 


1255.3 


1264.5 1265.4 


9.27 9.34 9.35 


561.0 565.1 


565.8 


Expt. c 


1.36075(3) 


1265.54(10) 


9.422(28) 


564, 604, 573 ± 42 m 


BeCl 


CVDZ 


1.8323 1.8269 1.8251 


817.4 


816.5 817.8 


5.12 4.89 4.89 


341.3 343.3 


344.7 


CVTZ 


1.8165 1.8117 1.8087 


833.1 


837.6 839.6 


4.87 4.90 4.93 


361.5 363.2 


364.8 


CVQZ 


1.8076 1.8018 1.7989 


839.3 


844.5 846.5 


4.91 4.94 4.94 


373.4 375.9 


377.4 


Expt. d 


1.7971 




846.7 


4.85(3) 


435, 385, 333, 384.8±9.2 m 


NaF 


CVDZ 


1.9869 1.9485 1.9478 


528.2 


527.2 527.7 


4.59 4.03 4.03 


449.4 455.4 


456.3 


CVTZ 


1.9966 1.9412 1.9402 


532.1 


523.9 524.4 


4.34 3.46 3.49 


462.1 466.1 


466.7 


CVQZ 


1.9967 1.9318 1.9311 


534.0 


530.4 530.6 


4.35 3.54 3.54 


469.0 473.2 


473.6 


CV5Z 


1.9916 1.9291 1.9284 


537.8 


533.8 534.0 


4.45 3.53 3.53 


473.1 475.4 


475.9 


Expt. e 


1.9259455(2) 


535.65805(21) 


3.57523(13) 


(514), 516 m 




NaCl 


CVDZ 


2.4278 2.4117 2.4099 


342.8 


345.9 346.8 


1.76 1.79 1.80 


383.9 384.4 


385.3 


CVTZ 


2.4133 2.3868 2.3845 


346.5 


354.2 354.4 


1.62 1.73 1.73 


396.2 394.2 


394.6 


CVQZ 


2.4064 2.3729 2.3705 


350.6 


360.2 360.7 


1.71 1.71 1.71 


407.5 404.8 


405.2 


Expt/ 


2.3607941(4) 


364.684163(391) 


1.776085(189) 


408, 409.3 ± 8 m 


MgF 


CVDZ 


1.7854 1.7782 1.7775 


697.9 


698.9 699.6 


4.19 4.18 4.18 


418.2 418.6 


419.7 


CVTZ 


1.7772 1.7627 1.7622 


701.0 


704.8 704.4 


4.00 3.96 3.94 


433.3 431.4 


431.9 


CVQZ 


1.7700 1.7541 1.7536 


708.4 


716.4 716.8 


4.07 4.20 4.21 


443.4 440.6 


441.0 


CV5Z 


1.7659 


711.7 




4.18 


447.1 




Expt.s 


1.7499371(1) 


720.14042(30) 


4.26018(16) 


458, 458.5 ± 5.0 
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TABLE VII: (continued) 





Val RIV All 


Val RIV All 


U) e X e 

Val RIV All 


D 

Val RIV All 


MgCl 


CVDZ 
CVTZ 
CVQZ 
Expt. h 


2.2582 2.2558 2.2541 
2.2251 2.2137 2.2110 
2.2193 2.2053 2.2025 
2.196(1) 


441.8 440.5 441.3 
457.6 461.9 463.2 
459.8 464.2 465.2 
466.0(8) 


1.85 1.79 1.79 
1.99 2.04 2.06 
2.01 2.04 2.05 
2.0(0) 


290.5 289.3 290.4 
310.5 307.1 308.1 
321.4 316.8 317.7 
317, 324.5 ± 2.1 m 


KF 


CVDZ 
CVTZ 
CVQZ 
CV5Z 
Expt. 1 


2.2436 2.2514 2.2506 
2.1993 2.1980 2.1972 
2.1868 2.1800 2.1796 

2.1714559(2) 


402.7 399.0 399.3 
414.7 415.6 415.7 
419.7 422.2 422.2 

426.261872(98) 


2.58 2.41 2.41 
2.46 2.24 2.26 
2.46 2.35 2.36 

2.449801(44) 


368.4 454.7 359.9 

332.5 477.6 370.1 
322.5 488.6 422.2 

319.6' 494.5* 435.6 ; 
489, 494.5 ± 2.5 m 


KC1 


CVDZ 
CVTZ 
CVQZ 
Expt/ 


2.9169 2.7709 2.7684 
2.9705 2.7045 2.7017 
2.9827 2.6829 2.6801 
2.666678(3) 


271.2 258.7 259.2 
301.9 270.4 270.8 
308.8 276.4 276.7 
280.07639(490) 


2.89 1.18 1.18 
3.31 1.15 1.15 
3.34 1.20 1.20 
1.31330(338) 


391.5 393.4 394.0 
400.4 407.8 408.3 
410.1 421.1 421.5 
419, 430.5 ± 8 m 


CaF 


CVDZ 

CV(D+3d)Z 

CVTZ 

CV(T+2d)Z 

CVQZ 

CV(Q+d)Z 

Expt.J 


2.1269 2.0406 2.0399 

1.9902 
2.1125 1.9889 1.9880 

1.9692 
1.9981 1.9597 1.9595 

1.9585 
1.9516403(1) 


405.4 536.4 536.5 

571.3 
399.2 562.3 563.0 

574.3 
537.6 583.6 583.5 

584.1 
588.644(2) 


30.03 2.62 2.62 

3.00 
38.60 2.75 2.78 
2.85 
5.37 2.87 2.87 
2.88 
2.91194(6) 


342.7 473.3 382.8 

508.5 
310.5 508.4 399.4 

532.3 
361.0 531.1 463.9 

531.6 
529, 524 ± 21 m 


CaCl 


CVDZ 

CV(D+3d)Z 

CVTZ 

CV(T+2d)Z 

CVQZ 

CV(Q+d)Z 

Expt. fc 


2.6598 2.5770 2.5752 

2.5049 
2.5969 2.4921 2.4891 

2.4704 
2.5402 2.4531 2.4501 

2.4519 
2.43674 


324.2 335.4 336.1 

348.2 
340.9 353.5 354.1 

357.1 
351.5 366.0 366.5 

366.4 
370.201 


1.22 1.16 1.15 

1.27 
1.27 1.27 1.27 

1.27 
1.34 1.35 1.35 

1.35 
1.3732 


24.8 352.8 353.7 
368.4 
370.0 382.0 382.8 

388.5 
396.7 404.1 404.9 
404.5 
395, 406 ± 9 m 



(a) (LiF) Hedderich, H.G., Engleman Jr., F.R., and Bernath, P.F., 1991, Can. J. Chem. 69, 1659. 

(b) (LiCl) Burkholder, J.B., Hammer, P.D., Howard, C.J., Maki, A.G., Thompson, G., and Chackerian Jr., C, 1987, J. Mol. 
Spectrosc. 124, 139. 

(c) (BcF) Tai, G. and Vcrma, R.D., 1995, J. Mol. Spectrosc. 173, 1. 

(d) (BeCl) Huber and Herzberg^. 

(c) (NaF) Muntianu, A., Guo, B., and Bernath, P.F., 1996, J. Mol. Spectrosc. 176, 274. 

(f) (NaCl, KC1) Sam, R.S., Dulick, M., Guo, B., Zhang, K.-Q., and Bernath, P.F., 1997, J. Mol. Spectrosc. 183, 360. 

(g) (MgF) Barber, B.E., Zhang, K.-Q., Guo, B., and Bernath, P.F., 1995, J. Mol. Spectrosc. 169, 583. 

(h) (MgCl) Rostas, J., Shafizadeh, N., Taieb, G., and Bourguignon, B., 1990, Chem. Phys. 142, 97. 

(i) (KF) Liu, M.-C, Muntianu, A., Zhang, K.-Q., Colarusso, P., and Bernath, P.F., 1996, J. Mol. Spectrosc. 180, 188. 
(j) (CaF) Kaledin, L.A., Bloch, J.C., McCarthy, M.C., and Field, R.W., 1999, J. Mol. Spectrosc. 197, 289. 

(k) (CaCl) Berg, L.-E., Klynning, L., Martin, H., Pereira, A2gmd Royen, P., 1981, Phys. Scr. 24, 23. 
(1) Single point energy calculation at experimental r e . 

(m) All except last values from Huber and Herzberg|^3l, last value from CRC Handbook]^] corrected to K with expt. io e x e 
and u> e . 



TABLE VIII: Computed and observed spectroscopic constants (kJ/mol, A, cm 1 as appropriate) 
for metal chalcogenides. 





Val RIV All 


Val RIV All 


LU e X e 

Val RIV All 


D 

Val RIV All 


BeO 


CVDZ 
CVTZ 
CVQZ 
CV5Z 
Expt. a 


1.3711 1.3674 1.3671 
1.3456 1.3414 1.3407 
1.3387 1.3334 1.3326 
1.3373 1.3317 1.3308 
1.3309 


1366.3 1371.1 1372.0 

1450.6 1457.1 1458.4 

1470.7 1483.5 1484.2 

1475.4 1488.2 1489.6 
1487.32 


11.47 11.22 11.18 
12.07 12.04 12.05 
12.07 12.11 12.12 
12.03 12.08 12.09 
11.830 


355.0 359.2 360.5 
408.5 412.7 413.6 
425.5 431.3 431.9 
429.9 435.8 436.6 
444, 431.0 ± 13. 4 d 


BeS 


CVDZ 
CVTZ 
CVQZ 
Expt. a 


1.7764 1.7704 1.7722 
1.7607 1.7568 1.7551 
1.7519 1.7467 1.7443 
1.7415(3) 


960.2 966.2 964.8 
977.5 982.6 983.7 
990.9 996.6 998.3 
997.94 


5.74 5.81 5.81 
5.82 5.83 5.86 
5.87 5.88 5.90 
6.137 


260.1 264.2 263.4 
290.4 293.0 293.8 

306.2 309.4 310.5 
367, 368 ± 59 d 


MgO 


CVDZ 
CVTZ 
CVQZ 
Expt. 6 


1.7859 1.7763 1.7753 
1.7662 1.7500 1.7491 
1.7595 1.7421 1.7413 
1.7481722(9) 


743.2 754.5 756.8 
772.2 794.1 795.8 
787.6 810.0 811.8 
785.2183(6) 


9.69 9.79 9.77 
9.21 9.04 8.94 
9.26 8.81 8.72 
5.1327(3) 


204.8 207.2 208.4 
234.0 236.5 237.0 
245.2 247.1 247.4 
(341), 359.7 ± 12. 6 d 


MgS 


CVDZ 
CVTZ 
CVQZ 
Expt. a 


2.1999 2.1952 2.1958 
2.1718 2.1600 2.1582 
2.1635 2.1491 2.1465 
2.1425 


499.9 498.6 498.4 
518.6 524.0 525.1 
524.4 530.3 531.5 
528.74 


2.63 2.51 2.50 
2.60 2.58 2.57 
2.62 2.60 2.59 
2.704 


162.5 163.3 163.8 

195.6 195.4 196.0 
208.5 207.4 208.0 

< 232, 231 d 



CaO 



CVDZ 


2.0360 2.0362 2.0531 


777.5 


779.9 


490.2 


2.17 


2.25 5.15 


143.5 283.0 194.2 


CV(D+d)Z 


2.0395 




557.1 






-10.94 


295.0 


CV(D+2d)Z 


1.9382 




477.8 






7.97 


308.0 


CV(D+3d)Z 


1.9095 




565.5 






11.78 


318.8 


CVTZ 


2.0206 1.9298 1.9284 


824.3 


523.1 


525.7 


31.47 


3.37 3.74 


63.0 332.5 230.3 


CV(T+d)Z 


1.8540 




661.2 






8.71 


360.5 


CV(T+2d)Z 


1.8472 




679.2 






7.88 


366.6 


CVQZ 


2.0658 1.8308 1.8307 


689.9 


718.0 


717.7 


5.23 


6.34 6.35 


66.8 390.7 331.4 


CV(Q+d)Z 


1.8281 




722.2 






6.13 


393.3 


CV(Q+2d)Z 


1.8274 




723.3 






6.11 


393.9 


DK-CVQZ 


1.8333 




709.6 






6.73 


385.0 


DK-CV(Q+d)Z 


1.8305 




714.1 






6.54 


387.5 


CV5Z 














406.4 e 


Expt. c 


1.8222315(4) 


732.01377(40) 


4.81268(74) 


> 459, 398.7 ± 16.7 d 



(a) (BeO, BeS, MgS) Huber and Herzbergl57l. 

(b) (MgO) Miirtz, P., Thiimmel, H., Pfelzer, O, and Urban, W., 1995, Mol. Phys. 86, 513. 

(c) (CaO) Focsa, O, Poclet, A., Pinchemel, B., Le Roy, R.J., and Bernath, P.F., 2000, J. Mol. Spectrosc. 203, 330. 

\~2 r~i 

(d) First value from Huber and Herzberg l57l . second from CRC Handbook 55] corrected to K with expt. u) e x e and oj e . 

(e) Single-point energy calculation at experimental r e . Results with the CV(5+d)Z and CV(5+2d)Z basis sets are 406.7 and 
406.6 kJ/mol, respectively. 



TABLE IX: Computed 2 D <- 2 S excitation energy (cm 1 ) in Ca + function of the basis 



D 



CVnZ 38065 22546 13965 12374 

CV(n+d)Z 28618 16539 12985 12898 

CV(n+2d)Z 21528 14920 12851 12885 

CV(n+3d)Z 19214 15155 
CV(Q+d)Z uncontracted 13240 
DK-CV(Q+d)Z uncontracted 14201 
Best estimate" 13846 



All calculations at the CAS-PT3 level l581 with a [Ne] core frozen. The experimental values are lsfll : 2 -D 3 y 2 13650.19, 
2 D 5/2 13710.88, spin-orbit averaged 13686.60 cm -1 . 

(a) CV(5+2d)Z plus relativistic correction estimated as the difference between DK-CV(Q+d)Z and CV(Q+d)Z. 
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